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Abstract 

A detailed analysis of dark matter event rates in minimal supergrav- 
ity models (MSGM) is given. It is shown analytically that the lightest 
neutralino the Z\ is the LSP over almost all of the parameter space, and 
hence the natural candidate for cold dark matter (CDM). The radiative 
breaking of SU{2) x U{1) constraints are shown to be crucial in deter- 
mining the expected event rates. Approximate analytic formulae are 
obtained to determine the gaugino-higgsino content of the Zi particle. 
Prom this one can deduce the behavior of the event rates as one varies 
the SUSY soft breaking parameters and tan f3. The constraint on the 
event rates due to the recently measured 5 — > s -|- 7 decay is calculated. 
It is seen that this data eliminates most of the parameter space where fi 
(the Higgs mixing parameter) and At (the t-quark cubic soft breaking 
parameter) have the same sign. Since the t-quark is close to its Landau 
pole, At is restricted to be mostly positive, and so most of the /x > 
part of the parameter space is eliminated. However, for /i < 0, one 
finds large regions of parameter space where the event rate is large and 
exceeds 0.01 events/kg da. The importance of proper treatment of the 
s-channel Z and Higgs poles in calculating the relic density is stressed. 
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The implications of the recent new experiments (SMC and EMS) on 
the quark polarizabihties are analysed and it is seen that uncertainties 
in these generally produce only small uncertainties in the event rates. 
A discussion is also given of the sensitivity of the expected event rates 
to changes in the allowed range of Z\ relic density. 
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1. Introduction 

The nature of the dark matter (DM) which makes up more than 90% 
of the matter of the universe is a particularly important issue as it may have 
a fundamental impact both on astronomy and particle physics. Dark matter 
has currently only been detected by its gravitational interactions, and thus 
may be composed of several constituents, e.g. baryonic dark matter (B), hot 
dark matter (HDM) and cold dark matter (CDM) (where "hot" and "cold" 
refer to whether the particle was relativistic or non-relativistic at the time of 
galaxy formation) . One may measure the amount of each species of dark mat- 
ter by the ratio Qj = Pi/ Pc where pi is the mass density of the ith constituent 
and Pc = SiJ^/SvrGiv is the critical mass density (H is the Hubble constant 
and Gjv is the Newtonian gravitational constant). Within the framework of 
the inflationary scenario one has = 1. The amount of baryonic matter 
is severely limited in the big bang cosmology by the observed abundancies of 
light elements i.e. Qb ~ 0.1. 

Rotation curves of stars imply a density of dark matter in our galaxy of 

Pdm = O.SGeV/cm^ (1) 

and this matter will be impinging on the Earth with velocity v^m — 320 km/s. 
The fact that microlensing find far more machos in the disk than in the halo 
of the Galaxy implies that at most 30% of the halo dark matter is machos [1]. 
Thus most of the halo of the Galaxy must be cold dark matter, and it is this 
dark matter that terrestial detectors can observe. 

A possible source of hot dark matter is massive neutrinos. In this paper we 
assume that the cold dark matter component is the lightest supersymmetry 
neutralino, the Zi particle. The anisotropy power spectrum (including the 
recent GOBE data) puts contraints on the relative amounts of HDM and GDM. 
A reasonable fit to the full spectrum gives flcDM '■ ^hdm = 2:1. Assuming 
fls — 0.1 one then estimates Q,^^ — 0.6. What is theoretically calculable is 
Qh'^ where h=(H/100 km/s Mpc). Current astronomical measurements yield 

0.4 ~ ~ 0.8 (2) 

i.e. two groups of measurements of h exist, one clustering at the lower bound 
and one at the upper bound. [The inflationary scenario (with zero cosmological 
constant) requires h ~ 0.5]. Thus we estimate 

0.1 < Qz^f^^ < 0.4 (3) 

and we will assume these bounds in the following. (Our results are not quali- 
tatively sensitive to the precise upper and lower limits of Eq. (3) and we will 
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discuss below what changes occur if one perturbs them.) We also note that 
it has recently been suggested [2] that if the value /i ~ 0.8 is correct, the age 
of the universe and other cosmological problems could be accounted for by a 
cosmological constant with Q ~ 0.6 — 0.8, and the remainder being CDM. This 
would also lead to Jl^i^^ being in the range of Eq. (3). 

Current dark matter detectors plan to obtain a sensitivity of i? ~ 0.1 
cvcnts/kg da. Future developments may improve this to R=0.01 events/kg 
da. We will thus limit our discussion here to the part of the parameter space 
where 

R ~ 0.01 events / kgda (4) 

since this sensitivity is what one may expect over the next 5-10 years. De- 
tection of the Zi depends on their scattering cross section by quarks in the 
nuclei of the detector. Thus a calculation of event rates depends on the two 
things: (i) that the relic density of Zi obey the bounds of Eq. (3), which limits 
the allowed SUSY parameter space, and (ii) a calculation of the Zi — q cross 
section. We consider these calculations in this paper within the framework of 
supergravity grand unification models [3] . While this model is not a complete 
theory it possesses a sufficient number of accomplishments to warrent using it 
as the dynamical framework. Thus it accounts naturally for grand unification 
at a scale Mq — 10^^ GeV implied by the LEP measurements of ai, 0L2 and 03 
at Mz\ it allows a natural breaking of supersymmetry (in the hidden sector) 
at the GUT scale (something that cannot be done in a phenomenologically 
acceptable way in the MSSM, and is yet to be demonstrated to occur in string 
theory); it can account for the suppression of FCNC interactions in a natu- 
ral way; in the minimal model (MSGM) it depends on only four additional 
parameters and one sign to describe all the masses and interactions of the 32 
new SUSY particles. (This may be compared with 110 new parameters that 
can occur in the MSSM.) 

The supergravity interactions of the MSGM produce four supersymme- 
try soft breaking terms at Mq scaled by (universal spin zero mass), mi/2 
(universal gaugino mass), and and (cubic and quadratic soft breaking 
constants). One of the remarkable features of this theory is that this sponta- 
neous breaking of supersymmetry at Mq generates, by radiative renormaliza- 
tion group (RG) corrections, the breaking of SU(2)xU(l) at the electroweak 
scale [4]: supersymmetry breaking produces SU(2)xU(l) breaking. We will 
see below that radiative breaking is a key element in the analysis of dark mat- 
ter event rates, and failure to include it loses much of the predictive power of 
the theory. 

While the MSGM possesses only four additional parameters to describe 
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SUSY phenomena, and this is far fewer than the 20 to 30 of the 110 possible 
new parameters that is commonly assumed in the MSSM, it still possesses a 
large parameter space. (Ideally, one would like to have four experiments to 
determine the four parameters, making all further predictions of the theory 
unique.) Recently, however, there have been two new pieces of data, the 
CLEO measurement of the 6 — > s + 7 branching ratio [5] and the CDF and 
DO measurements of the top quark mass [6]. We will see below that while 
large error flags still remain in these data, they greatly reduce the allowed 
parameter space. 

There has been considerable activity in the recent past to calculate ex- 
pected event rates for dark matter detectors [7-13]. However, Refs. [7-9] do 
not impose radiative breaking and thus can get abnormally high event rates 
(often by chosing the PC odd Higgs to be too light). The major part of the 
analysis of Ref. [10] is also done in this framework, and when radiative break- 
ing is introduced it is only for the special parameter choice Bo — Ao — rrio, 
and the entire parameter space is not scanned. Thus Ref. [10] predicts event 
rates that are too low. The analyses [7-11] also impose relic density constraints 
which leave out the full thermal averaging over the Z and h (light Higgs) s- 
channel poles. It is known that such omissions can generate serious errors 
in density calculations [14,15], and we will see below that it is important to 
treat s-channel poles correctly for about half the parameter space. Finally 
Refs. [7,8] do not include the heavy Higgs, H, in the event rate calculation, 
the importance of which was first pointed out in Ref. [16]. 

The plan of this paper is the following. In Sec. 2 we briefly review the 
ideas of radiative electroweak breaking and discuss the origin of the scahng 
relations between the masses of the light neutralinos and charginos and the 
gluino. In Sec. 3 we exhibit an approximate analytic formula for the gaugino 
and higgsino content of the Zi in the scaling regime. In Sec. 4 we discuss 
the relic density calculation, exhibit the importance of correct treatment of 
the s-channel resonances, and also discuss the event rate calculations. We also 
show that uncertainties in the nucleon spin content do not have any significant 
effect on event rates for all targets except for the lightest ones such as ^He 
and CaF2. These results are in contrast with a recent analysis [17] where the 
coherent part of the scattering was ignored and hence claimed a large effect. 
Sec. 5 is concerned with constraints on event rates and SUSY parameter space 
due to the 6 — > s -I- 7 decay and the top quark mass. We show there that there 
are sizable regions of the parameter space with R> 0.01 including this con- 
straint, and thus our results differ from those of Ref. [18] which concludes that 
with the 6 ^ s + 7 constraint, the event rate R < 0.01. We discuss in Sec. 6 
the effect of varying the endpoints of Eq. (3). Sec. 7 gives the conclusions. 
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The MSGM predicts that the Z\ is the hghtest supersymmetric particle over 
almost all the parameter space. The analytic analysis of this is given in the 
Appendix. 



2. Radiative Electroweak Breciking 

At the GUT scale the MSGM can be described by the superpotential 



where Wy is the cubic Yukawa couplings and W^'^ contains any quartic non- 
renormalizable couplings (which possibly lead to proton decay). The sponta- 
neous breaking of super symmetry leads to the soft supersymmetry breaking 
effective potential VgB and gaugino mass term L^^^^, 



where {za} are the scalar fields and A" the gaugino fields. Equations (5-7) 
arise after the supergravity interactions cause the spontaneous breaking of 
supersymmetry in the hidden sector, the GUT interactions cause the breaking 
of the GUT group G to SU(3)xSU(2)xU(l), all superheavy and hidden sector 
fields are integrated out, and non-renorahzable terms scaled by (1/Mp;) are 
neglected. The universality of the soft breaking parameters, mo, Ao, Bo and 
mi/2, is a consequence of the universality of gravitational couplings, plus the 
additional assumption that the hidden sector fields in the Kahler potential also 
couple universally to the physical fields. (This universality also quarantees 
suppression of FCNC interactions and thus is phenomenologically desirable) . 
The above results are generally insensitive to the nature of the GUT group 
provided the representations used to break G are not too large (so that GUT 
threshold corrections are not too big). 

The effective potential may be reduced to the electroweak scale by using 
the renormaliztion group (RG) equations. Minimizing the Higgs potential with 
respect to < i?i,2 > yields [4] 



(5) 



VsB = mlY^aZ^a + {AoWy + B0II0H1H2 + h.c) 
^mass = -mi/2X'^X°' 



(6) 
(7) 




Ej are loop corrections. The mjj. are given by 
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ml, = ml/^e{t) + A<,mi/2/(f) + mlh{t) - kit)Al (10) 

where the form factors e,f,g,h,k are defined in Ibanez et al [4], the gluino mass is 
mg — {o£3/aG)m-i/2 (with ac the GUT couphng constant), and t=ln(M^/(5^). 
One may show that solutions exist to Eqs. (8), i.e. that SU(2)xU(l) is 
spontaneously broken, if and only if at least one of the supersymmetry soft 
breaking interactions are non-zero. Thus it is the supergravity interactions at 
Mq that give rise to the breaking of SU(2)xU(l) at the electroweak scale. 

To obtain a qualitative picture of the implications of electroweak breaking, 
one notes from Eq. (8) that for tan? (5 >> 1 (i.e. tan(3 ~ 2 — 3) that 

~ -m]j^ - ^Ml - E2 (11) 

Thus for electroweak breaking to occur, it is necessary for to turn negative. 

(E2 is generally a small correction). The measured value of the top quark mass 
is mt = (176±8±10) GeV from CDF [6] and mj = (199t^5'±22) GcV from DO 
[6] while indirect determinations from LEP yield rrit = 16415*0+6-4 GeV [19]. 
These imply that the top is relatively close to its Landau pole. In this domain, 
for Q at the electroweak scale, h(t), e(t) and f(t) are negative and g(t) and k(t) 
are positive. Further, both | h \ and | e | are 0(1). (For a detailed discussion 
of the effects of the Landau pole see Rcf. [20].) Thus rnfj^ docs indeed turn 
negative allowing SU(2)xU(l) breaking to occur at the electroweak scale. 

In the following we will restrict nio and rrig to be less than 1 TeV to prevent 
an unreasonable amount of fine tuning. [Actually we will see that Eq. (4) 
actually implies rrig ~(650-700) GeV.]. Eq. (11) then implies 

/x^ >> M| (12) 

(which is satisfied if /i ~ (2 — 3)Mz) for almost the entire parameter space . 
Note also that the above discussion implies that /i^ is an increasing function 
of m? and m^. 

Eq. (12) is essentially the remnants of the gauge hierarchy problem in 
supersymmetry, i.e. from Eqs. (9-11), // is scaled by rrio and rrig and we allow 
the latter to go as high as 1 TeV. It has been previously shown [21] that Eq. 
(12) leads to a set of scaling relations between the light neutralinos, chargino 
and gluino, 

1 1 

2mz^ - rriz, = m^^ - (3 " 1)^-9 (^3) 
as well as the additional relations 

"^Z3,4 - "^1^2 "^^i^ ruH = mH± = mA» ruh. (14) 
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Also one finds 

63Gey <mh^ mCeV (15) 

where the lower bound on mh is the current LEP hmit. Eqs. (12-14) are thus 

a direct consequence of radiative electroweak breaking, and we will see below 
that they play a dominant role in determining relic densities and dark matter 
detection rates. 



3. Composition of Zi 



The Zi is generally a mixture of gauginos W^jB and Higgsinos Hi,H2. 
We write 

Zi = niWs + + risHi + n^H2 (16) 

The expansion coefficients are determined by diagonalizing the neutralino 
mass matrix. In the (W3, B,Hi,H2) basis this reads [3] 



where = {ai/asjrrig, a 

—MzsinOwcosP and d = M, 



/ m2 


a 


^ \ 


rhi 


c 


d 


a c 







\ h d 


-1^ 


J 



(17) 



= Mzcosdwcos(3, b — —Mzcos9wsinf3, c — 
zsinOwsinjS. The Zi is the lowest mass eigen- 



vector of M^, which may be determined in general numerically. However, as 
discussed in Sec. 2, radiative breaking implies >> Mz- One may generate 
an approximate analytic form by treating the two blocks proportional to Mz 
perturbatively. To second order perturbation theory one finds 



ni 



n2 



n3 = 



^4 



IM, 



2 II {1 — ml/ iM^) m2 — rfii 

sin^O 



sin29w 



sin2p + 



^_ 1M| 



2/^2 (1 - mf///2)2 

— ^ 9 , ^ sinOwsinP 

//I — mil n'^ 



w 



mi 
1^ ' 



1 + —ctn(3 



M. 



1 



fjL 1 — fhyiJ,'- 



:sin9wCOs(3 



mi ^ 
1 H tanp 



mi 


(18) 


. 




(19) 








(20) 




(21) 



Eqs. (18-21) differ from the numerical computer results by amounts Srii ~ 0.03 
over almost the entire parameter space (and are generally a good deal better). 
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Thus one may use them to understand the nature of the solutions. One sees 
first that the Zx is mostly bino since n2 deviates from unity by second order 
effects, 0{M\I Usually, n2 > 0.95 and often larger. However, 713,714 and 
n\ are first order, 0{Mzl This allows the Higgsino components of the Z\ 
to become considerable e.g. ris 0.2. The coherrent part of the Zi-nucleus 
scattering in dark matter detectors, depend upon the interference between the 
gaugino and Higgsino components of the Zy [7]. Thus such terms can become 
quite large, and as will be discussed in Sec. 4, this means that the coherrent 
scattering almost always dominates the incoherrent (spin dependent) scatter- 
ing. Thus there is a large difference between a Z\ whose bino amphtude is 
0.95 and one which is 100% bino, and one cannot approximate the former as 
being pure bino. We also note that the above results are a direct consequence 
of the radiative breaking conditions, and would not in general hold without 
them. 



4. Relic Density and Event Rate Analysis 

The primordial Zi, created at the time of the Big Bang, can annihilate in 
the early universe. The main diagrams for this are shown in Fig. 1 [22]. At 
high temperature, the Z\ is in thermal equilibrium with its decay products. 
However, when the annihilation rate falls below the expansion rate of the uni- 
verse, freeze out occurs at temperature Tj. The Z^ are then disconnected from 
the background and continue to annihilate. In the simplest approximation 
(accurate to about ± 25 %) the current relic density is given by 




^ 2.5 X 10-1^ ( ^ ) ( ^ ) (22) 
where 

J{xj) = dx < av >; Xf ~ kTf/m^;^ (23) 

Jo 

Here x=kT/m^^, Nj is the effective number of degrees of freedom at freeze 
out, {T^_ /T.y)^ is the reheating factor, and the current microwave radiation 
temperature (in °K). In Eq. (23), a is the annihilation cross section, v is the 
relative velocity and < > means thermal average. 

Freeze out generally occurs when the Zi are non-relativisitic, i.e. ai Xf — 
kTf/m^_^ ~ ^. Thus the thermal average may be taken over a Boltzmann 
distribution: 



< av >-- 



/ dvv^av exp[—v^/Ax\/ / dvv^ exp[—v^ / Ax\ (24) 
Jo Jo 
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The non-relativistic nature of the annihilation process has lead, in the past, 
to making a non-relativistic expansion of av i.e. av = a -\- h{v'^ / c^) H — •; after 
which the thermal average becomes trivial to take. This is a good approxima- 
tion for the t-channel pole diagrams of Fig. 1, and for the s-channel diagrams 
when 2m^^ is not in the vicinity of the h or Z poles. However, when 2m^^ 
is near the s-channel poles, the non-relativistic approximation fails [14] and 
can produce errors as large a factor of 100 [15] due to the fact that crv is a 
rapidly varying function in this region. Further, the thermal averaging smears 
the region where this effect can occur (characteristically over a region pa 10 
GeV in m^^ or f« 50 GeV in nig). For this situation one can first perform the 
integral of Eq. (24) analytically, and then calculate Eq. (23) numerically. 

In the following, we will investigate the parameter space defined by the 
following bounds: 

lOOGel/ <mo< ITeV; IbOGeV < m-g < ITeV (25) 

- 6 < At/nio < 6; 2 < tan(3 < 20 (26) 

The lower bound on m^ is the current Tevatron bound and the upper bounds 
on m,,, m^ are to prevent excessive fine tuning of parameters. (Similarly, tan/3 
~ 30 represents a fine tuning of the Higgs VEV ratio.) The range on Aj 
covers the allowed parameter space (when all other experimental constraints 
are imposed). One may estimate the region of this parameter space where the 
above discussion of s-channel poles is important. The scaling relations, Eqs. 
(13,14), allow us to choose m^ as the independent variable. One sees from 
these relations, that, varying over the field parameter space, one is generally 
near an s-channel resonance (h or Z pole) when m^ ~ 450 GeV, and hence for 
this region one must treat the s-channel terms accurately. This is borne out 
by detailed numerical calculations which show that significant errors occur 
in the calculation of ^Zi^-^ with nig < 450 GeV when the non-relativistic 
approximation to < av >is made, while the approximation is generally good 
for m^ > 450 GeV. Since the constraint on detection rate, R> 0.01 events/kg 
da, requires mg~ (650-700) GeV, we see that a correct treatment of the s- 
channel poles is important for a large fraction of the total parameter space. 

Detection of dark matter impinging on the Earth depends upon the Zi- 
quark scattering cross section for the quarks in the nuclei of the detector. The 
basic diagrams are shown in Fig. 2. This process has been studied by a number 
of authors [24] , and can be represented by the following effective Lagrangian: 

Leff = {Xil^l''Xi)[(llM<iPL + B,PR)q] + {xiXi){(lC,m,q) (27) 
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Here xi is the Zi field, q is the quark field and Pr,l — (l±7^)/2. The coef- 
ficients Aq and arise from the Z t-channel pole and the q s-channel pole, 
while Cq comes from the h° and H° t-channel poles and the q s-channel pole. 
(We follow the notation of Ellis and Flores [7] where explicit formulae for Aq, 
Bq, Cq are given for the q, h and Z pole diagrams.) 

The first term of Eq. (27) gives rise to spin dependent (incoherrent) scatter- 
ing while the second term gives rise to spin independent (coherrent) scattering. 
Summing over all the quarks in the nucleus, the latter term then produces an 
additional factor of the nuclear mass. Mat. The event rate for a detector then 
takes the form [24] 



R — [Rsi + Rsd] 



Q.SGeVcm-^ 



320km/s 



events 
kg da 



(28) 



where p^.^ is the local mass density of ZijV2^ the incident Zi velocity and the 
spin independent (SI) and spin dependent (SD) rates have the form 



R 



SI 



MM + m 



Zi 



Asi? 



Rsd 



IGrriy Mat „ , , , , 
.^^^X^J{J + l)\AsD 



Mjv + rriz^ 



(29) 
(30) 



where J is the nuclear spin and A is defined by<A'^|X^5'i|-/V>=A<A^| j| 
>. We note that for large Mat, R5/ ~ Mn while Rsd ~ I/^a^j showing that 
the heavy nuclei arc best for detectors sensitive to spin independent scattering. 

Since the squarks are generally heavy over most of the parameter space, 
A57 is dominated by its Higgs contributions which has the general form [26] 



tHiggs 
^SI 



9l 


41 


( cosau 1 
sinfi 




sinaj] 
sinfi 


1] 


U 


— quark 


4Mw 


mi 


—sinaH 
L cos/3 J 


mjj 


cosaH 
y cosp J 




d 


— quark 



(31) 



where Fh — {rii — n2tan6w){niCosaH +nzsinaH) and Fh = {rii — n2tan9w) 
{n4sinaH — n^cosaH) and an is the rotation angle needed to diagonalizc the 
h-H° mass matrix. Thus the SI scattering arises from interference between the 
gaugino and Higgsino parts of the Zi i.e. from Eqs. (18-21) from the n2 x 
terms of F^^h for most of the parameter space. In general (including the loop 
corrections to the Higgs mass matrix [27]) one finds that tanan ~ 0.1. Hence, 
for most of the parameter space, the h contribution to d-quark scattering is 
suppressed by a factor of tan^a/f relative to the H, and this can overcome the 
fact that mn "^/mh ^ >> 1. In fact we find that the H contribution varies 
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from 1/10 to 10 times the h contribution as one varies over the full parameter 
space, and it is essential to keep both neutral CP even Higgs bosons in the 
analysis [16]. 

In contrast to the above, the contributions to Asd from the Z pole depends 
on nl — n\ which is small by Eqs. (18-21). Thus for heavy nuclei detectors 
one always has 

Rsi » RsD (32) 

and even for the hght CaF2 detector which has a large \^J{J + 1) value, the 
spin independent scattering dominates over most of the parameter space. 
^SD depends on the spin content of the nucleons defined by 

< n I ^7^? I n >= 2sf^)Ag (33) 

where s^^-^ is the spin 4-vector of nucleon n, and Aq measures the part of 
the nucleon spin carried by quark q. It has been suggested [17] that the 
differences between older data [28] and more recent data [29] determinations 
of Aq (particularly As) could produce uncertainties in the value of R leading 
to errors as large as a factor of 30. This could indeed be the case if the Zi 
was pure gaugino (i.e. n2 = 1, na = n4 = rii = 0) for then R5/ would vanish. 
However, as we saw above, a significant interference between the gaugino and 
Higgsino parts of the Zi exists (i.e. the Zi amplitude has a gaugino amplitude 
of only about 0.95) leading instead to Eq. (32) for heavy nuclei, minimizing 
the effect. Fig. 3 shows the ratio between the predicted value of the total R 
using the new and old data for Aq for a Ge detector. The difference is less 
than 10% over the entire parameter space. While larger errors can exist for a 
CaF2 detector, where R5D is large, even here the difference is less than ± 30% 
over 92% of the parameter space. 

A more serious uncertainty exists in R5/ due to a lack of knowledge of the 
s-quark content of a nucleon defined by < n | russs \ n >= fsM^- Estimates 
of fs [30] have about a 50% uncertainty leading to a ± (30-50) % uncertainty 
in R57. This will not, however, change the qualitative nature of the results 
given below. 

One may now understand analytically the dependence of the event rate on 
the different SUSY parameters. From Eq.(ll), we see that //^ is an increasing 
function of -mf^^ hence by Eq. (10) an increasing function of m^ and mp. 
As shown in the discussion following Eq. (31), the major contribution to A5/ 
(which dominates the total event rate R) is proportional to (71-2^3) [or (n2n4)] 
which by Eqs. (18-21) are the leading [0{M.z/ lA terms. These decrease with 
increasing 11. Thus one expects R to be a decreasing function of m^ and m^. 
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This indeed was what was seen in the detailed computer calculations of Ref. 
[12]. Further, for radiative breaking, which implies tan > 1, Eq. (31) also 
shows that the d-quark amplitude is an increasing function of tanj3. Again the 
computer calculations of Ref. [12] show this rapid rise of R with tan {3. 

The calculation of the detector event rates now proceeds as follows. One 
calculates ^zj^^ selects that part of the parameter space of Eqs. (25,26) 
that satisfies the constraint of Eq. (3) as well as the LEP and Tevatron bounds 
on SUSY masses. One then calculates the event rate R for this allowed part of 
the parameter space. Figs. 4 and 5 show the maximum and minimum event 
rates for // < and // > as a function of At/rrio for Ge and CaF2 detectors 
as one lets the remaining parameters vary over the allowed parameter space. 
(The relevant parts of these graphs are for the regions where R„ia.T > 0.01.) 
As can be seen from the previous discussion, the largest event rates occur for 
the largest allowed values of tan P and for the smallest values of m^. (The 
sharp peaks and dips in the Rmax curves arise from the fact that a small value 
of m^ implies by Eq. (13) a small m^|.^. If this parameter point also satisfies 
the LEP cut that m^^^ > 45 GeV it is allowed and gives rise to a large event 
rate. Otherwise it is excluded.) The Ge detector is generally considerably 
more sensitive then the CaF2 detector since its nuclear mass is considerably 
larger increasing the Rsi contribution [as seen from Eqs. (29, 30)]. 

Figs. 4 and 5 also show that the parameter space is bounded in At with 
most of the allowed aregion occurring forA^ > 0. This phenomenon is due 
mainly to the fact that the t quark mass is large and hence close to its Landau 
pole. As one approaches the Landau pole, the light stop mass, m^^ obeys [20] 

where 

AR = At-0.613m-g; = 0A6A[{m{f - {mty]/M^ (35) 

m{ = 197sinP is the fixed point mass, and m?^ (NP) is a relatively smooth 
non-pole contribution. For A|. sufficiently large, the ti becomes tachyonic, 
eliminating such parameter points from the parameter space. One expects 
then a lower bound on A^ for A^ < 0, and similarly an upper bound for A^ > 
(as we are requiring m^ < 1 TeV), with the domain of positive Aj being 
larger than for the negative At since cancelation in Ar between the A^ and m^ 
term can occur in the former case. One finds, in fact, including correctly the 
non-pole part in Eq. (34), that for mj ~ 170 GeV 

- 1.5 ~ At/ruo ~ 5.5 (36) 
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Thus the high mass of the t quark ehminates a large amount of the SUSY 
paramter space. 



5. Constraints Prom 6 — > s + 7 Decay 

The 6 — > s + 7 decay is a sensitive test for new physics, since it is a flavor 
changing neutral current (FCNC) processs. Thus the Standard Model and 
new physics loops enter at the same order. This is explicitly exhibited in Fig. 
6 where the basic diagrams for the decay at scale ~ Mw are given for SUSY 
models. The W-t loop is common with the Standard Model, while the other 
loops are the additional supersymmetric contributions. The measured CLEO 
branching ratio for S — > + 7 is [5] 

BR{B Xs-f)) = (2.32 ± 0.5 ± 0.29 ± 0.32) x 10"^ (37) 

or combining all errors in quadrature one has BR(B^ -^s7) — (2.32 ± 0.66) x 
10~^. This result represents an additional limitation on the allowed SUSY 
parameter space, and we discuss in this section the effect the CLEO data has 
on dark matter detection event rates. 

In the spectator approximation, the B meson decay can be related to the 
b quark decay. It is convenient, to define the quantity R as 



BR(B XcePe) T(b^c + e + Pe) 

(A CKM and (m^)^ factor, which have large errors, cancels out in R.). Here 
BR(B^ XcePe) = (10.7± 0.5)%. The diagrams of Fig. 6 can be described by 
an effective Hamiltonian [31] 

Heff = Vt,VZ^Cr{Mw)Q7 (39) 

where Q7 = {e/2AT^'^)mi,SL(T'^'^hRF^j,. Here F^,^ is the electromagnetic field 
strength and rri}, is the b-quark mass. One must use the renormalization group 
equations to go from scale = My/ to /i ~ where the decay occurs. To 
leading order (LO) in QCD corrections, R is then calculated to be [31] 

^H^I^Jm"^'"*-'!^ '''' 

where I{z) = 1 — 82;^ + 82;^ — — is a phase space factor (z=mc/mb) 

and 
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C'/f(m,) = rj^sCr(Mw) + -(?7^ - V^')C8(Mw) + C^iM^) (41) 

Here Qs = {g3/167i^)mhSRa'^^T^bLG'^^, where and G^^ are the gluon gener- 
ators and field strengths and C2 represents operator mixing with the 4-quark 
operators. 

The QCD corrections arc large for this process and the next to leading order 
correction (NLO) are needed to obtain an accurate theoretical prediction. At 
present, however, not all the NLO terms have been calculated. Thus the 
theoretical analysis has an estimated error of about ± 30% with a Standard 
Model (SM) prediction of BR[B X^^] ^ (2.9 ± 0.8) x 10"^ for m^ = 174 
GeV [32]. The H~-t SUSY diagram of Fig. 6 adds constructively to the SM 
amplitude while the W — i loop may enter constructively or destructively with 
the SM model amplitude. Since the central value of the CLEO data of Eq. 
(37) already lies about 1 std. below the central value of the SM results, the 
current data cannot tolerate a large amount of constructive interference with 
SUSY amplitudes. Thus, in spite of the large errors in the current data, the 
fe — s + 7 decay produces a significant constraint on the SUSY parameter 
space. 

As seen from Eqs. (14), radiative breaking generally implies mH± is large, 
suppressing its contribution to the 6 — > s + 7 decay. The i —W diagram can 
become large, however, when m^^ and m^^ are small. The t {massY matrix 
reads 



mt{At + iictnl3) 



mt{At + ^ctnP) 



(42) 



where expressions for m|^, are given in [20]. The light t eigenvalue is 

(43) 



2 -'- / 2 I 2 \ 



^{ml - m}^f + ml{At + ^Jictn(5f 



Thus m?^ can become small if Aj and ji can have the same sign and this effect 
is enhanced due to the fact that m^ is so large [33] . One expects therefore that 
when Kt and have the same sign the theoretical prediction for BR(b— > s-l-7) 
will become large if also m^^ is small, while the SUSY effect on BR (b^ s + 7) 
will be small for A^ and ^ having opposite sign (or mj^,^ becoming large). This 
can be seen explicitly to be the case from detailed computer calculations given 
in [34]. 

Eq. (36) shows that the major part of the parameter space has At > 0, 
and so the b— > s + 7 decay is expected to infiuence mostly the // > branch. 
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One may quantify this by requiring that the theoretical rate for BR(b— > s + 7) 
be within the 95 % CL bounds of the experimental value. One finds then that 
for > one is restricted to the region 



i.e. about 40% of the parameter space of Eq. (36) is eliminated by the b— s+7 
data! In addition, sections of the parameter space with small nig are eliminated 
for At < 0, /X < and for < At < 0.5, > 0. 

The effects of the above restrictions are shown in Fig. 7 (/x < 0) and Fig. 
8 (/X > 0). In Fig. 7 one sees that event rates arc reduced in the small domain 
of At < 0, (as the low m^^ part of the parameter space is eliminated), but 
not significantly modified over the remainder of the parameter space where 
At > 0, since here // and At have opposite signs. Fig. 8 shows that the only 
remaining part of the parameter space is the narrow band of At < 0.5, and all 
of the parameter space with At > 0.5 (where fi and At have the same sign) is 
eliminated. Thus the major effect of the b^ s + 7 data is to eliminate regions 
of parameter space, rather than reduce expected event rates. As can be seen 
from Figs. 7 and 8, there still remain regions of parameter space with large 
event rates. (This is to be contrasted with the results stated in [18].) 

6. Varying The Bounds On f^^^/i^ 

The analysis considered above was done within the framework of the bounds 
of Eq. (3) on VL^^h?. We consider now the effect of varying these bounds. We 
first note that the Zi annihilation cross section in the early universe arising 
from the diagrams of Fig. 1 is a decreasing function of m^^. Alternately, from 
the scaling relations Eq. (13), one may say then that ^Zi^"^ increases as m^^ 
increases. On the other hand, we saw in Sec. 4 that the event rate R is a de- 
creasing function of m^. We are interested in this paper in the region R >0.01 
events/kg da which will be accessible experimentally in the forseeable future. 
This bound on R then puts an upper bound on m^. The largest allowed value 
of m^ occurs for the largest value of tan/3 and smallest value of m^ (which by 
Eqs. (25,26) wc arc taking as tan,'?< 20, nio > 100 GcV), since R increases 
with tan/3 and decreases with nio. One finds then for the parameter space 
defined by Eqs. (25,26) that for < one has[35] 



For II > 0, lUg can rise to (700-750) GeV. However, this occurs for At > 0, and 
as discussed in Sec. 5, this region of the parameter space is eliminated by the 



1.0 < At/nio < 0.5; > 



(44) 



m~g ~ 650^61^; for R > 0.01, // < 



(45) 
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b— > s + 7 decay data, and so Eq. (45) represents the true upper bound on 
for dark matter that can be detected with current designs of detectors. 

A detailed inspection of the fuU parameter space shows that the hmitation 
R< 0.01 imphes then that the early universe annihilation channels Z\Z\ 
hh, Zh are closed. The channel Z^Zi -^WW,TL is almost always closed, and 
the small regions in parameter space which allow this annihilation are very 
close to threshold and hence highly suppressed. Hence, one need not consider 
the vector meson channels in the analysis of ^zj^"^ given in Sec. 4. 

The behavior oiVL^^h? and R as a function of m^ also shows that the bound 
R> 0.01 implies that mg < 650 GeV. This is exhibited in Fig. 9a. One sees, as 
shown in Fig. 9b, that ^Zi^^ <0.3 when R> 0.01 at these maximum values of 
rrig. Thus in this domain, the results obtained above are not sensitive to the 
precise upper bound in Eq. (3). If one were to lower the upper bound below 
Q^^/i,2 = 0.3, then the upper bound on m^ is further reduced. This is shown 
in Fig. 10, where it is seen that if f^^i^^ < 0-2) then m^ < 400 GeV, which 
would make the gluino accessible to the proposed high luminosity upgrade of 
the Tevatron [i.e. the TeV(33)]. We note that in both the inflationary scenario 
with a cold-hot dark matter mixture or the scenario with cold dark matter and 
a cosmological constant, low values of flzi^^ preferred. This follows from 
the fact that in the former case one needs a small value of h (i.e. h ~ 0.5 
so that the age of universe be consistent with the estimated age of globular 
clusters) and in the latter case because Jl^^ is small (i.e. Q^i — 0-2 — 0.4 since 
the majority of the matter of the universe is in the cosmological constant). 

We now turn to the question of sensitivity of results to the lower bound on 
Qz^h'^- As discussed above, low ^Zi^^ arises when m^ is small, which is also 
the domain of parameter space where R can be large. Also, as discussed in 
Sec. 4, the peaks and dips of Rmax in Figs. 4 and 5 arise from whether or not 
the LEP bound m^y^ > 45 GeV can be satisfied. Figs. 11 and 12 exhibit the 
effect on the maximum event rates for Pb and CaF2 detectors when the bound 
^Zi^^ > 0.10 is raised to Jl^^/i^ > 0.15. One sees that the sharp peaks get 
chpped off, but otherwise the results are qualitatively unchanged. There still 
remains, however, a sizable amount parameter space where R exceeds 0.01. 

7. Conclusions 

We have examined in this paper the direct detection possibility of Zi cold 
dark matter within the framework of the minimal supergravity model (MSGM) 
for the parameter space defined by Eqs. (25,26), and for relic densities ^z^h^ 
in the range given by Eq. (3). [Eq. (3) encompasses the range one would 
expect from infiationary cosmology for either the cold-hot dark matter senario 
or the cold dark matter plus cosmological constant possibility.] One further 
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limits the parameter space so that experimental bounds on SUSY masses are 
obeyed. Two new pieces of data, the t quark mass and the b— > s + 7 branching 
ratio, have greatly constrained the SUSY parameter space. Thus the fact that 
uit is large (i.e. close to its quasi infra-red fixed point) limits the domain of At 
to be mostly positive, while the experimental b— > s + 7 decay rate limits the 
parameter space to be mostly where and /i have the opposite sign. Thus 
the majority (though not all) of the allowed parameter space is in the region 
where > and fj, < 0. 

Physical quantities in the MSGM depend on four SUSY parameters, m^,, 
TXig, At, tan/3, and the sign of /i. Thus a general quantity has a complex 
behavior as one varies over the full parameter space. Radiative breaking, 
however, plays a central role in MSGM predictions and allow one to understand 
analytically the qualitative behavior of the event rate for the detection of Zi 
dark matter. Approximate analytic expressions were obtained in Sec. 3 for 
the content of the Zi showing that the Zi was mostly bino, but with a non- 
negligible amount of higgsino. One can see from this that the spin independent 
(cohcrrcnt) contribution generally dominates R (and hence the most sensitve 
detectors are those with the heaviest nuclei) and that R decreased with m^, 
and mg and increases with tan/?. 

Current dark matter detectors hope to obtain a sensitivity of R> 0.01 
events/kg da in the forseeable future. As seen here, such a sensitivity will al- 
low the exploration of a sizable amount of the SUSY parameter space, though 
there will still remain large sections that fall below this bound. The domain R 
>0.01 correpsonds to m^ ~ 650 GeV and is insensitive to the choice of upper 

bound on ^Zi^^ '^^ ^'i- (3) provided f^^i^^ ^ 0-3) ^^'^ oiily mildly sensitive to 
the lower bound of Eq. (3). Thus dark matter could only be expected to be 
seen by current detector designs if the gluino is not too heavy. 

Acknowledgements 

This research was supported in part by NSF grant numbers PHY-9411543, 
PHY-19306906 and PHY94-07194. 



18 



Appendix 

In N=l supergravity there are three possible candidates for the LSP. These 
are the Zi, the v and the cr. We discuss here the region of parameter space 
where the Zi is the LSP. The basic requirement for this is then that 

m^j,>mz^] ■mo>mz^ {A.l) 

In addition one has the LEP constraint that 

rriz^ ~ 2QGeV {A.2) 

The requirement that the be heavier than Zi is an experimental constraint 
as a charged LSP would have already been discovered. We will see that in fact 
this condition is obeyed for almost the entire parameter space. The region 
where the D is the LSP will also be seen to be very small. 

The cr and masses can be expressed in terms of the basic MSGM pa- 
rameters as [4] 

mj^ = ml + aalfiml - stn^ewM^cos2p {A.3) 

m\ + |M|cos2/3 {AA) 
where the form factors are defined as 

(2 + Ai)/(l + At)2 (A5) 

where t = 2in{MG/Q) and A = aoif, 1, - 3) are the U{1)XSU{2)XSU{3) 
/3-f unctions with aa — ao/^Tr. In the following we use ac = -^^Mg = 2 x 
10^^GeV,Q = Mz. One finds then /i ^ 38.1, /s = 99.0, /g = 772.0. In 
order to analytically illustrate the results, we will perform the remainder of 
the calculation in the scaling limit of Eq. (13,14). (A more accurate numerical 
calculation give results close to the analytic ones.) One can then relate mi to 

the Zi mass by 

=7711 = (q;i/q;g)?7t.i (^-6) 
or mi = 2.45m^^. Eqs. (A. 3, A. 4) then become 

m?^ ^ml + 0.912m|^ - sin'^ewM^cos2p {A.7) 

ml^ml + 3.19m| + |M|cos2/? {A.8) 
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The requirement that the cr not be the LSP then becomes an upper bound 



on m^^: 

m^^ < 11.4 [ml + siv?ewMl{-cos2(3)] (A.9) 

Similarly, the condition that the v is heavier than the Zi becomes a lower 
bound on m^^ 

> 0.456 [-ml + iM|(-cos2/3)] (A.IO) 

Eqs. (A. 9, A. 10, A. 2) are three inequalities constraining the parameters m^^, 
mo and tan/3. In addition we have the fine tuning constraints of m^,, m^ < 
1 TeV. 

In the scaling limit then Eq. (A. 9) implies for m^ = 1 TeV 

mo > 41.9, for tanP = 1; m<, > 0, for tanp > 4.7 (A.ll) 

with corresponding smaller lower bounds on m,, for smaller values of m^. Thus 
for almost the entire parameter space the MSGM predicts that the is not 
the LSP. 

Turning to the u constraint, we see that the r.h.s. of Eq. (A. 10) falls below 
the LEP bound Eq. (A. 2) (and hence becomes irrelevant) when 

ml > \Ml{-cos2(3) - (20)70.456 {A.12) 

and hence 

mo > for tanP < 1.24; > 57.3^6^ for tanp » 1 {A.13) 

Alternately for all tan/3, Eq.(A.lO) is satisfied for all mo if m^^ > 43.6 GeV 
(or uig > 308 GeV) . Thus the z> is predicted to be heavier than the Zi for all 
but a very small portion of the parameter space. 
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Figure Captions 

Fig. 1 Z\ annihilation diagrams for annihilation in the early universe. 

Fig. 2 Zi-quark scattering diagrams for a terrestial dark matter de- 
tector. 

Fig. 3 R(New)/R(01d) vs \x for a Ge detector. "New" data is Ref. 
[29] and "Old" data Ref. [28]. 

Fig. 4 Maximum and minimum event rates for Pb (solid) and CaF2 
(dashed) detectors as a function of Aj for ^l <0 and m^ = 168 GeV. The 
other parameters are varied over the range of Eqs. (25,26). (From Ref. 
[12].) 

Fig. 5 Same as Fig. 4 for > 0. 

Fig. 6 Diagrams contribution to b^ s + 7 decay at W mass scale. 

Fig. 7 Maximum and minimum event rates for Pb detector vs Aj/mo 
without b— > 5 + 7 constraint (solid) and with b^ s + 7 constraint (dot- 
dashed) at 95% CL for /i < 0. Other parameters are varied over the 
range of Eqs. (25,26) with mt=168 GeV. (The minimum event rates are 
unaffected by the b— > s -|- 7 constraint.) 

Fig. 8 Same as Fig. 7 for // > 0. The parameter space with the 
6 — > s -|- 7 constraint terminates for Kt/rrio > 0.5. 

Fig. 9 Maximum value of m^ (Fig. 9a) and maximum value of ^z-J^^ 
at this mg(Fig. 9b) as a function of tan/3 for < for domain R> 0.01 
events/kg da for Pb detector as other parameters are varied over range 
of Eqs. (25,26) with m^ = 168 GeV. The curves from bottom to top are 
for kt/mo = 0, 2.0, 4.0. 

Fig. 10 Maximum value of m^ as a function of the upper bound on 
Jl^i^^ for // < for tanl3 — 6, kt/nio — 0.5. (Results are insensitive to 
the values of tanP and Aj.) 

Fig. 11 Maximum event rates for Pb detector as a function of At/rrio 
for /i < for fi^i^^ > 0.10 (solid) and fi^^/i^ > 0.15 (dot-dashed). Other 
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parameters are varied over the range of Eq. (25,26) with uit —168 GeV. 
(The b— > s + 7 decay constraint is not here imposed.) 

Fig. 12 Maximum event rates for CaF2 detector as a function of 
At /mo for /X < for Q^i^^ > 0-10 (dashed) and fi^i^^ > 0-15 (sohd). 
Other parameters are varied over the range of Eqs. (25,26) with m^ = 
168 GeV. (The b— > s + 7 decay constraint is not here imposed.) 
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